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Abstract 



In this paper we consider some flexible discrete sampling of a discrete scale in- 
variant process {X(t),t 6 R + } with scale I > 1. By this method we plan to have 
q samples at arbitrary points So, Si, . . . , s g _i in interval [I, I) and proceed our sam- 
pling in the intervals [l n ,l n+1 ) at points i n So, l n &i, • . . , / n s g _i, n £ Z. Thus we 
have a discrete time scale invariant (DT-SI) process and introduce an embedded 
DT-SI process as W(nq + k) = X(l n s k ), q G N, k = 0, . . . , q - 1. We also con- 
sider V(n) = (V°(n), . . . , where V k (n) = W(nq + k), as an embedded 
q-dimensional discrete time self-similar (DT-SS) process. By introducing quasi Lam- 
perti transformation, we find spectral representation of such process and its spectral 
density matrix is given. Finally by imposing wide sense Markov property for W(-) 
and V(-), we show that the spectral density matrix of V(-) and spectral density 
function of W(-) can be characterized by {Rj(l), Rj(0), j — 0, . . . ,q — 1} where 
R j (k) = E[W(j + k)W(j)}. 
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similar processes. 

1 Introduction 

The concept of stationarity and self-similarity are used as a fundamental property to 
handle many natural phenomena. Lamperti transformation defines a one to one corre- 
spondence between stationary and self-similar processes. Discrete scale invariance (DSI) 
process can be denned as the Lamperti transform of periodically correlated (PC) process. 
Many critical systems, like statistical physics, textures in geophysics, network traffic and 
image processing can be interpreted by these processes pQ. Fourier transform is known 
as a suited representation for stationarity, but not for self-similarity. A harmonic like 
representation of self-similar process is introduced by using Mellin transform [1]. 
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A process which is Markov and self-similar, is called self-similar Markov process. 
These processes are involved in various parts of probability theory, such as branching 
processes and fragmentation theory [2]. 

Current authors considered DSI processes in the wide sense with some scale I > 1. 
They proposed to have some fixed number of samples, say T, in each scale at points 
a k , k G Z where / = a T , T € N. By such sampling they provided a discrete time scale 
invariant process in the wide sense and found a closed formula for its covariance function 
[UJ. In this paper we consider X(-) as DSI process with scale I > 1, and sampling at 
arbitrary points 1 ^ So < s± . . . < s g _i < I in the interval [1,1). We also take our samples 
at points l n s$, l n si, . . . , Ps g _i, n € Z in the intervals [l n , l n+1 ). Then we introduce some 
discrete time embedded scale invariant (DT-ESI) process W(nq + k) = X(l n Sk), q € N, 
k = 0, . . . , q— 1 and corresponding multi-dimensional discrete time embedded self-similar 
(DT-ESS) process as V(n) = (V°(n), . . . ,F 9_1 (n)) where V k (n) = W(nq + k). We 
investigate properties of these processes when they are Markov in the wide sense. 

This paper is organized as follows. In section 2, we present a review of multi- 
dimensional stationary, periodically correlated, self-similar and discrete scale invariant 
processes. Then we define discrete time self-similar (DT-SS) and scale invariant (DT-SI) 
processes. We also introduce quasi Lamperti transformation in this section. Section 3 is 
devoted to the structure of the multi-dimensional DT-SS process resulting from the above 
method of sampling. We define DT-ESI process and corresponding multi-dimensional 
DT-ESS process and characterize the spectral density matrix of it in this section. Finally 
covariance function and spectral density matrix of the discrete time embedded scale in- 
variant Markov (DT-ESIM) processes and corresponding multi-dimensional discrete time 
embedded self-similar Markov (DT-ESSM) are obtained in section 4. 

2 Theoretical framework 

This section is organized in tree subsections. First we review the structure of the co- 
variance function and spectral distribution matrix of multi-dimensional stationary pro- 
cesses. We present definitions of DT-SS, DT-SI, wide sense self-similar and scale invariant 
processes in subsection 2.2. We define quasi Lamperti transformation and present its 
properties which provide a one to one correspondence between DT-SS and discrete time 
stationary processes and also between DT-SI and DT-PC processes. 

2.1 Stationary and multi-dimensional stationary processes 
Definition 2.1 A process {Y(t),t G R} is said to be stationary, if for any t, r 6 R 

{Y (t + T)}±{Y(t)} (2.1) 

where = is the equality of all finite- dimensional distributions. 
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If (2.1) holds for some r G R, the process is said to be periodically correlated. The 
smallest of such r is called period of the process. 

By Rozanov [8], if Y(t) = {Y k (t)} k =1 ...n be an n-dimensional stationary process, 
then 

Y(t) = [ e iXt (j)(dX) (2.2) 



is its spectral representation, where (f) = {<p k }k=l ... n an d '■Pk is the random spectral 
measure associated with the fcth component Y k of the n-dimensional process Y. Let 



B kr (r) = E[Y*(t + t)Yr{t)], k, r = 1, . . . , n 

and B{t) = [-Bfcr (■?")] fc,r=l,...,n be the correlation matrix of Y. The components of the 
correlation matrix of the process Y can be represented as 

B kr (r) = J e iXr F kr (d\), k, r = 1, . . . , n (2.3) 



where for any Borel set A, F kr (A) = E[(p k (A)(p r (A)] are the complex valued set functions 
which are cr-additive and have bounded variation. For any k,r = 1, . . . ,n, if the sets 
A and A' do not intersect, E{(p k (A)ip r (A')] = 0. For any interval A = (Ai,A2) when 
F kr ({Xi}) = F kr ({\2}) = the following relation holds 

oo 

F fcr (A) = — / Y, B k r(T)e~ iXT d\ (2.4) 

J A T=-00 

1 1 p -i^2T _ p -iXiT 

= —B kr (0)[X 2 - Ax] +lim- V B kr {r) : 

2ir T^oo 2tt ^— ' —it 

0<\r\^T 

in the discrete parameter case, and 

1 f a e~ iXiT — e~ iXlT 
F kr (A) = lim — / : B kr (r)dT 

a^co 2TT J_ a —IT 

in the continuous parameter case. 

2.2 Discrete time scale invariant processes 

Definition 2.2 A process {X(t),t G R- + } is said to be self-similar of index H > 0, if 
for any A > 

{X- H X(Xt)} £ {X(t)}. (2.5) 

The process is said to be DSI of index H and scaling factor Ao > or (H,Ao)-DSI ; if 
(2.5) holds for A = Aq- 
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As an intuition, self-similarity refers to an invariance with respect to any dilation 
factor. However, this may be a too strong requirement for capturing in situations that 
scaling properties are only observed for some preferred dilation factors. 

Definition 2.3 A process {X(k),k € T} is called discrete time self-similar (DT-SS) 
process with parameter space T, where T is any subset of countable distinct points of 
positive real numbers, if for any k\ , k 2 £ T 

{X(k 2 )} I (h)« { x( kl )}. (2.6) 
fci 

The process X(-) is called discrete time scale invariance (DT-SI) with scale I > and 
parameter space T, if for any ki,k 2 = lk\ € T, (2.6) holds. 

Remark 2.1 If the process {X(t),t € R + } is DSI with scale I = a T for fixed T £ N 
and a > I, then by sampling of the process at points a k , k € Z, we have X(-) as a DT-SI 
process with parameter space T = {a k , k € Z} and scale I = at T . If we consider sampling 
of X(-) at points a nT+k , n € Z,/or fixed k = 0, 1, . . . , T — 1, then X(-) is a DT-SS process 
with parameter space T = {a nT+k ,n € Z}. 

Yazici et.al. [9J introduced wide sense self-similar processes as the following definition, 
which can be obtained by applying the Lamperti transformation Ch to the class of wide- 
sense stationary processes. This class encompasses all strictly self-similar processes with 
finite variance, including Gaussian processes such as fractional Brownian motion but no 
other alpha-stable processes. 

Definition 2.4 A random process {X(t),t € R + } is said to be wide sense self-similar 
with index H, for some H > if the following properties are satisfied for each c > 0, 

Mi>*2 > o 

(i) E[X 2 (t)] < oo, 

(ii) E[X(ct)\ = c H E[X(t)\, 

(Hi) E[X(ctx)X{ct 2 )} = c 2H E[X{tx)X{t 2 )]. 

This process is called wide sense DSI of index H and scaling factor cq > 0, if the above 
conditions hold for some c = c$. 

Definition 2.5 A random process {X(k), k € X 1 } is called DT-SS in the wide sense with 
index H > and with parameter space T, where T is any subset of distinct countable 
points of positive real numbers, if for all k,k\ € T and all c > 0, where ck,ck\ € T : 

(*) E[X 2 (k)] < oo, 

(ii) E[X(ck)} = c H E[X(k)}, 

(Hi) E[X(ck)X(cki)] =c 2H E[X(k)X(k 1 )}. 
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If the above conditions hold for some fixed c = cq, then the process is called DT-SI 
in the wide sense with scale cq. 

Remark 2.2 Let {X(t),t G R + } in Remark 2.1 be DSI in the wide sense, with the same 
scale I = a T . Then X(-) with parameter space T = {a k ,k G Z} for a > 1 is DT-SI 
in the wide sense and X(-) with parameter space T = {a nT+k ,n G Z} for fixed TgN, 
a > 1 is DT-SS for k = 0, . . . , T — 1 in the wide sense. 

Through this paper we are dealt with wide sense self-similar and wide sense scale 
invariant process, and for simplicity we omit the term "in the wide sense" hereafter. 

2.3 Quasi Lamperti transformation 

We introduce the quasi Lamperti transformation and its properties by followings. 

Definition 2.6 The quasi Lamperti transform with positive index H and a > 1, denoted 
by £>H,a operates on a random process {Y(t),t G R} as 

C H ,aY(t) = t H Y(\og a t) (2.7) 

and the corresponding inverse quasi Lamperti transform Cj^ a on process {X(t),t G R + } 
acts as 

C- H ) a X{t) = a~ tH X{a t ). (2.8) 

Corollary 2.1 If {Y(t),t G R} is stationary process, its quasi Lamperti transform 
{£H, a Y(t),t G R + } is self-similar. Conversely if {X(t),t G R + } is self-similar pro- 
cess, its inverse quasi Lamperti transform {£^ a X(t),t G R} is stationary. 

Corollary 2.2 // {X(t),t G R+} is (H,a T )-BSI then C H ) a X{t) = Y(t) is PC with 
period T > 0. Conversely if {Y(t),t G R} is PC with period T then Cu,aY{t) = X(t) is 
(H, a T )-DSI. 

Remark 2.3 If X(-) is a DT-SS process with parameter space f = {l n ,n G Z}, then its 
stationary counterpart Y(-) has parameter space T = {nT, n G Z} 

X(l n ) = C H , a Y{l n ) = l nH Y{\og a a nT ) = a nTH Y(nT). 

Also it is clear by the following relation that if X(-) is a DT-SI process with scale I = a T , 
T G N and parameter space T = {a n ,n G Z}, then Y(-) is a discrete time periodically 
correlated (DT-PC) process with period T and parameter space T = {n, n G Z} 

Y{n) = C H ] a X{n) = a- nH X{a n ). 
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3 Structure of the process 



In this section we define a multi-dimensional DT-SS process in the wide sense. We also 
introduce a new method for sampling of a DSI process with scale I > 1, which pro- 
vide sampling at arbitrary points in the interval [1,1) and at multiple l n of such points 
in the intervals [l n ,l n+1 ), n 6 N. We introduce DT-ESI process corresponding to the 
multi-dimensional DT-ESS process. Finally in Theorem 3.1 we find harmonic like rep- 
resentation and spectral density matrix of the multi-dimensional DT-ESS process. 



Definition 3.1 The process U{t) = (U°(t), U 1 ^), . . . , U q 1 (t)) with parameter space 
f = {l n , n G Z} ; / = a T , a > 1 and T £ N is a q- dimensional discrete time self similar 
process in the wide sense, where 

(a) {U 3 (-)} for all j = 0, 1, • • • , q — 1 is DT-SS process with parameter space 

fi = {l n ,n£ Z}. 

(b) For every n, r G Z, j, k = 0, 1, ■ ■ ■ , q — 1 

Cov([P(P +r ),[/ fe (P)) = l 2nH Cov(U j (l T ),U k (l)). 

Our method of sampling is to provide enough flexibility to choose arbitrary sample 
points of a discrete time scale invariant process. So, one could decide to have T partitions 
in each scale interval I n = [l n ,l n+1 ), n G Z of a continuous time DSI process X(-) with 
scale / > 1 and find a by / = or . Then our partitions in scale interval are 

[a nT , a nT+1 ), [a nT+1 ,a nT+2 ), [a nT+T ~\ a^ T ). 

So we consider to have nj, samples in partition [a nT+k , a nT+k+1 ) at points 

where 1 ^ s kl < s k2 < ■ ■ ■ < s knk < a, k = 0, . . . , T — 1 and q = Y2i=o n «- Now we can 
state the following remark. 

Remark 3.1 Let U k (l n ) = X(l n s u ) in Definition 3.1, where s u = a k s kx in which 

Yli=-i n i ^ u < Yli=-i n i> n -i = an d u = x + Yli=o n i> x = 1, Thus 
X(l n s u ) foru = 0,...,q-l is a DT-SS process and U(l n ) = [X(l n s ), 
is a q-dimensional DT-SS process. 



By such method of sampling at discrete points we provide a (/-dimensional DT-ESS 
process V(n) as 

V(n) = (V°{n),V 1 (n),...,V q - 1 (n)), n€ Z 
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where q = Yli=o n i an< ^ 

V u {n) := X(a nT s u ) (3.1) 
Yli=-i n-i^u < Ya=-i n i-> n -i = 0, s u = a k s kx and x = u- YaZI n i, u = 0, . . . , q - 1. 

Remark 3.2 Corresponding to the q- dimensional DT-ESS process V(n) there exist a 
DT-ESI process W(k) with scale I = a T as 

W(k) := V u (n) = X(a nT s u ) k G Z (3.2) 

where u = k — q[^], n = [^] and k = nq + u, since by (3.1) and (3.2) 

W(k + q) = X(o> +1 ) T s u ) = a TH X{a nT s u ) = I h W{k). 

By the following theorem, the spectral density matrix of the (/-dimensional DT-ESS 
process and harmonic like representation of each column is obtained. 

Theorem 3.1 Let X(-) be a DSI process with scale I = a T and 1 ^ so < Si < . . . < 

s 9 _i < a T , then V(n) = (V°(n), ... , V q ~ l (n)), where V u {n) = X(a nT s u ), n G Z and 
u = 0, . . . , q — 1 is a q- dimensional DT-ESS process and 

(i) The harmonic like representation ofV u {n) is 

V u (n) = (a nT s u ) H 1^ j" n d<\> u (oS) (3-3) 
Jo 

where (j) u {w) is an orthogonal spectral measure, that is E[d(j) u (oj)d(f) u (oj')\ = 0, u, v = 
0, . . . , q — 1 when oj / oo' . 

(ii) The corresponding spectral density matrix ofV(n) is g H (oj) = [g^,(w)]«,z/=o,...,cj-i> 
where 

r= — oo 

r G N and Qu yU { T ) * s ^ e covariance function ofV u (r) and V u (0). 

Proof of (i): Remark 2.3 implies that 

V u (n) = X(a nT s u ) = £ H , a Y(a nT s u ) = (a nT s u ) H V u (n) 

where n u (n) = Y(nT +log a s u ). Thus V u (n) for every u = 0, l,...,q— lisa DT-ESS pro- 
cess in n, where its discrete time stationary counterpart rj u (n) for fixed u = 0, 1, . . . , q — 1 
has spectral representation 7? M (n) = J Q 27r e tujn d(f) u (oj). 
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Proof of (ii): The covariance matrix of V(n) is denoted by Q H (n, r) = [Q^,(n, T)]„ jl/= o,..., g -i 
where 

Q£„(n,r) = £[F> + r)F»] = E[X(a^ T s u )X{a nT s v )] 

By the scale invariant property of the process X(-) we have that 

Q% v (n, t) = a 2nTH E[X(a^ T s u )X(s u )] = a 2nTH Q^ u (r) (3.5) 

where Q« v {t) = Q*,(0,r) = £[^(^"(0)], then by (3.3) 

f'2ir flit 

Q£„(t) = E[(cJ T s u ) H {s v ) H / e^dcp u ^) / dM^ 7 )] 

Jo Jo 

= a^ TH (s u s u ) H e^eZGf » (3.6) 

JO 



where £ , [<i(/> u (£j)d(/> I/ (a;')] = dG^ v {uj) when uj = uj' and is when w 7^ uj' . 

On the other hand, by the definition of rf (n) in the proof of part (z) 

Q^(r) = £[X(a TT s u )X(s,,)] = fi^/ta^s^if/fs,)] 
= (q^s^I^tT + log a s u )Y(log a s„)] 
= (a^s u s v ) H E[t(rWm = (a^ s u s u ) H B UyU (r). 

Then by (3.6) 

,-2-k 

B u , u (r)= e^dG^iu), u,u = 0, ...,q- 1 

JO 

Now by (2.3) and (2.4) for u, v = 0, . . . , q — 1 we have that 

1 r 00 

G£„(A) = -/ £ ^(rJe-^dA. 

r=-oo 

By substituting B Uj1/ (t) = {a TT s u s v )~ H u (t) , the elements of the spectral distribution 
function, G^ u (-) has the following representation 

= / ^ a"™^"^* (r)dA. (3.7) 



Let A = (u,u) + du>], then the elements of the spectral density matrix, g^ u (uo) are 







T— — OO 



Thus we get to the assertion of part (ii) of the theorem. □ 



4 Multi-dimensional DT-ESSM process 

Using our method of sampling in section 3, we find the covariance function of the DT-ESI 
process W(-), which is defined in (3.2) and its corresponding multi-dimensional DT-ESS 
process V(-), defined in (3.1) for the case that they are Markov in the wide sense as well, 
which we call them DT-ESIM and DT-ESSM respectively in subsection 4.1. We find the 
spectral density matrix of these processes in subsection 4.2. 

4.1 Covariance function of DT-ESIM 

Here we characterize the covariance function of the DT-ESIM process {W(k),k € Z} 
in Theorem 4.1 and the covariance function of the associated q-dimensional DT-ESSM 
process in Theorem 4.2. 

Theorem 4.1 Let {W(k), k € Z}, defined in (3.2), be DT-ESI and Markov in the wide 
sense DT-ESIM, with scale or . Then for r G W = {0, 1, . . .}, k = nq+v, k+t = mq+u, 
u, v = 0, . . . , q — 1 and n, m 6 Z, the covariance function 



R k (t) := E[W{k + t)W(k)} = E[X(a ml s u )X(a nl s„)] 



(4.1) 



can be characterized as 



Rn(tq + s) = [f(q - l)]*/> + s - 1)[/(k - l)]" 1 ^) 
R K (-tq + s)= a~ 2tqH R K+S {(t -l)q + q-s) 



(4.2) 



where 1 ^ sq < Si < . . . < s y _i < a , t G Z, s = 0, . . . , q 



1 



r r 



/»=n/w=n^( i )/^(°)' 



r G Z 



(4.3) 



and /(-l) = 1. 



9 



Before proceeding to the proof of the theorem we present the main property of covariance 
function of the wide sense Markov process. 

Let {X(n),n G Z} be a second order process of centered random variables, E[X(n)] = 
and E^A^n)) 2 ] < oo, n € Z. Following Doob [3], the real valued second order process 
X(-) is Markov in the wide sense if 

R(n\,n2) = G( min(ni, n2))H ( max(ni, 112)) (4-4) 

where R{n\,n2) := E[X(ni)X(ri2)} is the covariance function of X(-) and G and H are 
defined uniquely up to a constant multiple and the ratio G/H is a positive nondecreasing 
function. 



Proof of the theorem: As {W(k),k G Z} is DT-ESI with scale a T , this theorem 
fully characterize the covariance function of the DT-ESIM process. From the Markov 
property (4.4), R k (t) defined in (4.1), satisfies 

R K ( T ) = G{a nT s u )H(a mT s u ), r G Z, a > 1 (4.5) 
By substituting r = in the above relation we have m = n, then 



G{a nI s„ 



Therefore 



Thus 



= fro^ '- T 6 z (4 - 6) 

where H(k) = H(a nT s u ) and H(k + r) = H(a mT s u ). Therefore 

R.(r) = ^^^.(0). (4-7) 



So 



For r = 1, we have 



By the recursive relation, it follows that 

R K {l)R K ^(] 
R K (0) R K -i(0) "' Rq(0) 



H <^^l€^-ti H (»)- H (°)ri/o) 



j=Q 
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and 

K-l 

H(«) = h(o) n M 

3=0 

where f(j) = Rj(l)/Rj(0). By the assumptions n = [^], v = K — q[^] we have k = nq + v, 
then 

nq+v— 1 

H(n<z + i/) = H(0) jj /(j). (4.8) 

As mentioned in Remark 3.2, {W(k),k G Z} is DT-ESI with scale /, then 

_ _ E[W(k + g + 1)W(k + g)] 

+ i? K+g (0) E[W(K + q)W(K + q)] 

_ a 2TH E[W(k + = = 

a 2 ™£[W( K )If( fi )] i? K (0) 

Hence by (4.8) 

H(ng + ,/) = H(0) [ J] /(*)] " II " > 1 

By the definition of / in (4.3) 



H(nq + v) = H(0)[/(g - !)]"/> - 1). (4.9) 



By a similar method one can verify that 

H(-nq + u) = H(0)[/(g - l)]- n /> - 1). 



Let r = tq + s in (4.7), t € W and s = 0, 1, . . . , q — 1, then it follows from (4.9) that 

V ; H(/c) kV ' H(0)/(k-1) 

For r = — + s we have that 

+ s) = J B[X(a-' 9+K+s )X(a K )] = a - 2tqH E[X{a K+s )X{a tq+K )} 
= a~ 2t " H R K+S (tq -s) = a- 2t « H R K+s ((t -l)q + q- s).D 



Now we can use this theorem to prove the next result as follows, for g-dimensional 
DT-ESSM process. 
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Theorem 4.2 Let {W(k),k G Z} be a DT-ESIM process, and {V(n),n G Z} be its 
associated q- dimensional DT-ESSM process with covariance matrix Q H (n,r) which is 
defined by (3.5). Then 

Q H (n, t) = a 2nTH [f(q - 1)} T CR, r G Z (4.10) 

where /(•) is defined in (4.3) and the matrices C and R are given by C = [C Uji ,] U)iy= o : .., i? -i, 
where C UjU = f(u — l)[/(f — and R is a diagonal matrix with diagonal elements 

R u (0), v = 0, 1, . . . , q — 1, which is defined in (4.1). 

Proof: As W(-) is DT-ESI with scale I, (3.2) and (3.5) indicate that Q^ u (n,r) = 
a 2nTH qH ^^ T y ]\f ow ^he assumption k = nq + v and k + t = mq + u where m, n G Z, 
r G W, we have r = (m — n)q + u — v and therefore 

R k (t) = R nq+U ((m - n)q + u - v) = E[W(mq + u)W(nq + v)] 

= E[X(a mT s u )X(a nT s u )]. 

Hence 

Qu,Ar) = E[X(a TT s u )X( Su )} = R v (rq + u-v) (4.11) 

and by the Markov property of W(-) from (4.2) we have 

R u (rq + u- v) = [f(q - l)] T }\u - l)[/> - l)]" 1 ^) 
for u, v = 0, . . . , q — 1. Let C UjU = f(u — 1)[f{v — l)] -1 , so 

Q^(r) = [f(q-1)YC U;U R V (0). (4-12) 

Thus we can represent the elements of the covariance matrix of g-dimensional DT-ESSM 
process as 

Q^(n,r)=a 2nTH [f(q-l)VC u ,M0)n 
4.2 Spectral representation of the process 

The spectral density matrix of the g-dimensional DT-ESSM process is characterized by 
the following lemma which is proved in [TJ. 

Lemma 4.1 The spectral density matrix g H (oj) = [g^ v (oj)]u,u=o,...,q-l of the q- dimensional 
DT-ESSM process V(n) is specified by 

/(n-l)fi.(O) fjy - 1)^(0) 

/> - 1)(1 - e-^a- HT f(q - 1)) f(u - l)(l - e-^a HT f- x (q - 1)) 

where -Rfc(O) is the variance ofW(k) and /(•) is defined by (4.3). 
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Proof: By applying (3.4) and (4.12), the spectral density matrix of the process {V(n),n £ 
Z} which is denoted by g H {ui) = [g^ v {u))\ u ,v=o,...,q-\ can be written as 



H i, \ { s uSy) I" —THT —iuiT/~iH / \ 

f u ,z/H = — ^ — [2^ a e v«,^l r ) 

1 r=0 
+ £ «- THr e— Q^(r)] = <&>) + <&>) 

r= — oo 

where 

Ci(-) = 1!li | L - E a ~ THTe ~ tu)T ^ - wfa - !)[/> - ^r 1 ^ ) 

T=0 

._ {*u*u)- H Ku ~ l)ik(O) 



2tt/(i/-1) 



T = 



By Remark 3.2, the scale invariant property of W{k) and the assumption, that at least 
one of the Con:\W(j)W(j + 1)] be smaller than one, we have that \f(q — 1)| < a TH for 
j = 0, . . . , q — 1. Thus 

\e-^a- TH f(q - 1)| = |a- TJ 7(g - 1)| < 1, 
and (4.13) for r € W is convergent. By the equality 

QuA-r) = E[X( a - TT s u )X(s u )} = a~ 2TTH E[X(a TT s u )X(s u )} = a~ 2rTH Q U)U ( T )> 
convergence of g^ v 2 (w) follows by a similar method. Therefore 

»u,i/v y 2vr L /(^ _ l) ^ v w ;y 

i4(0)/>-l) 



_^5>-™ e -/ ((? -l )r 



/(« ^ .=i 

(s u s,)- H [■ «,(0)/(n - 1) + _R„(0)/> - l)a- TH e^f(q - 1) 



2vr L/(„ - 1) (1 - a~ TH e-^f{q - 1)) f(u - l)(l - a-™e^f{q - 1)) 
so we arrive at the assertion of the lemma. □ 



Example 4.1 Xei 

oo 
n=l 
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where B(-) is the standard Brownian motion, /(•) indicator function, H > and A > 1. 
We call this process Simple Brownian Motion. We showed in JW that {X{t),t £ R + } is 
DSI and Markov with Hurst index H and scale A. By sampling of this process at points 
a s u , n € W, where 1 < So < s%, ■ ■ ■ , s 9 _i < a T , and by assuming A = a T , 

W(k) : = X(a nT s u ), 

is a D T-ESIM process, and V(n) = (V°(n), . . . , V q -i(n)) where V u (n) = W(k) is the 
associated q-dimensional DT-ESSM process where u = k — q[^], n = [£]. By (4-1) we 
have that Rf(0) = Rf(l) = a 2TH 'sj for j = 0, • • • ,q-2 and R^l) = a TH ' R^O) = 
a 3TH 's q ^, where H' = H — ~. So R u (0) = a 2TH 's u , R u (0) = a 2TH 's u . Also (4.3) 
implies that f(u — 1) = f{y — 1) = 1, f(q — 1) = a TH . Thus By Lemma 4-1, the spectral 
density matrix of V{n) is 

1 - e~ iul a- T l 2 1 - e- iuJ a T / 2 _' 
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